# Full text of "Magnetic topology and current channels in plasmas with toroidal current density inversions"

## See other formats

CO o CN Or < 43 c*2 o Oh' > On CN O o en s Magnetic topology and current channels in plasmas with toroidal current density inversions D. Circfl and I. L. CaldasQ Departamento de Fisica Aplicada, Universidade de Sao Paulo, 05508-090, Sao Paulo, Brazil. The equilibrium magnetic field inside axisymmetric plasmas with inversions on the toroidal cur- rent density is studied. Structurally stable non-nested magnetic surfaces are considered. For any inversion in the internal current density the magnetic families define several positive current chan- nels about a central negative one. A general expression relating the positive and negative currents is derived in terms of a topological anisotropy parameter. Next, an analytical local solution for the poloidal magnetic flux is derived and shown compatible with current hollow magnetic pitch measurements shown in the literature. Finally, the analytical solution exhibits non-nested magnetic families with positive anisotropy, indicating that the current inside the positive channels have at least twice the magnitude of the central one. PACS numbers: 52.55.-s, 52.30.Cv, 41.20.Gz In tokamak devices a toroidal magnetic field confine the orbits of the charged particles inside a chamber. A toroidal current flowing within the plasma produces a complementary magnetic field preventing particle drift looses. Non-inductive drive mechanisms help to sustain this current for long pulses with small or negative induc- tive drive and during a slow transition from positive to negative toroidal current [1,] . A relevant question arising from this situations is that of the structure of the mag- netic field if the toroidal current density becomes neg- ative in some region of the plasma. In the last decade, the achievement of quasi-steady-state alternating current scenarios [2| and the observation of stiff structures with nearly zero magnetic pitch angles in a finite region about the plasma center [3, |j] has attracted attention to the problem of current reversal equilibrium configurations (CRECs). In this work we study the equilibrium topology of the magnetic field subjected to azimuthal (toroidal) current density inversions in an axisymmetric plasma. In the fol- lowing pages it is shown that the equilibrium topology is composed of non-nested families of nested magnetic sur- faces, where each magnetic family defines a current chan- nel inside the plasma. The relation between the currents in the channels is studied in terms of topological quan- tities. In a more quantitative approach, an analytical solution of the equilibrium problem about a region of in- terest provides the topology of the magnetic surfaces and the relevant control parameters as well as their bifurca- tion values controlling the transition between different equilibrium configurations. This is done without speci- fying further plasma profiles or arbitrary functions. The obtained solution agrees with several published equilib- ria [5|-l8| while providing a simple understanding of the control parameters and their relation with the current in the different channels. For axisymmetric systems the equilibrium magnetic field may be written B = W x V^ + FV(j>, (1) where i/j(R,z) = RA^R, z) and F(R,z) = RB^{R,z) are proportional to the poloidal magnetic flux and cur- rent respectively. B$(R,z) and A$(R,z) are the az- imuthal components of the magnetic field and vector potential and (i?, 4>, z) the usual cylindrical coordinates. The term "Vip x V(f> corresponds to the poloidal magnetic field B pi i.e. the projection of the magnetic field over a plane <j> = const. It is clear that B ■ V-0 = 0, so that the magnetic field lines remain attached to the level sets of ip(R, z) called magnetic surfaces. For axisymmetric magnetic confinement devices the main source of poloidal field is the plasma current, con- sequently, the level sets of i/){R, z) form families of nested tori inside the plasma. In equilibrium, the toroidal cur- rent flowing inside a magnetic surface labeled by ip is given by the Ampere's law as fM)It(ip)=f B-dl = (j>-f Vij)xdL = ±j> \Vip\dl, Jr^p j^tp -^r^, (2) where T^ is a magnetic circuit resulting from the inter- section of the magnetic surface ip with an arbitrary plane <f> = const. For a simple magnetic circuit the product V"0 x dl always points in the same direction, the current is negative when V?/> is inwards and positive otherwise. The existence of current density reversals suggest that some magnetic surface ip® contains a vanishing current; from (J5J this requires |V-0| = along the continuous magnetic circuit Tq = T^, . This requires that for any co- ordinate system {u, v}, the equation d u ip(u, v) = leads to the same relation between u and v that the equation d v ip(u,v) = 0. In addition, such relation must satisfy ip(u,v) — tpa- This kind of degeneracy is possible for one-dimensional problems, but leads to structural insta- bility in general two-dimensional equilibria [7|, i.e. any variation in the poloidal field, will destroy the topology of the surface ip = tpo. A more feasible situation, i.e. structurally stable, re- quires that VV> vanishes at isolated points inside the plasma. Some of these points must be saddle, each one introducing four branches where the condition ift(R, z) — ip c is satisfied, with tp c the value of the poloidal flux on the critical point. This guarantees that the poloidal field reverses over the surface ip(R, z) = ip c after vanishing at each saddle. The CRECs topologies may be classified depending on how the separatrix r c = Tv, , connects two branches of the same hyperbolic point [9] . If T c connects two oppo- site branches of a saddle the resulting circuit is simple but non-unique, i.e. there is a second simple circuit T' c connecting the remaining two branches; this leads to an even number of axisymmetric islands (Fig. QJ,). On the other hand, when T c connects two non-opposite branches the circuit is non-simple but unique and leads to an odd number of islands (Fig. [TJj). Separatrixes delimit sev- eral families of nested magnetic surfaces acting as current channels inside the plasma. These axisymmetric chan- nels arc also called non-nested surfaces in the literature. In Fig. Q] the poloidal field direction reveals the sign of the toroidal current inside each channel. In summary, any internal toroidal current density reversal in a two- dimensional magnetic equilibrium leads to a chain of ax- isymmetric islands with positive current about a central structure with negative current. simple circuits -{Ti,^} (Fig. [2]d) it can be shown that FIG. 1. Odd and even systems of axisymmetric islands. Cor- responding diamonds are connected by the separatrix F c . The arrows indicate the direction of B p and " ■" or " x" in the mag- netic axes indicate the current direction for each channel. As the poloidal field only vanishes at isolated points and never reverses in a regular magnetic surface, its line integral on a magnetic circuit can never be zero. From this, the current enclosed by a magnetic circuit just out- side the positive channels is positive, in consequence, the total current in the positive channels must exceed the current of the negative one, not only equate its magni- tude. Each positive channel contains an elliptic point and is bounded by a separatrix and two hyperbolic points [lOj ■ We can always build a curve 7 parallel to V-0 that passes through all the hyperbolic and elliptic points defining the positive current channels (Fig. [2^,). The circuit 7 encloses a vanishing current, given that B p -dl = in all its points. By decomposing the critic magnetic circuit r c into two / B p di = J2 r x = h , i B p di = Y^r 2 = i 2 , (3) where I\ 2 is the current flowing through the region be- tween 7 and T12 on the i'th positive current channel. The relative difference between the half currents on the i'th channel is 7?; = (I| - I\)/I\ with /* = /» + I| the total current inside that channel. The current flowing inside all the positive channels is I + = I\ + I2 and the one flowing in the central channel is /_ = —I\. Using the previous relations we have F 1-Tj' (4) with 77 = YliVi^i/Y^i^i- Equation (fj| relates the total current in the positive channels with the negative cur- rent in the central one through the parameter 77, that measures the anisotropy of the positive current channels. For small anisotropies we have i+ ~ 2J_ and in simple cases where the islands are created in regions with mono- tonic variation of the current density we expect r\ > 0, leading to 1+ > 2|J_|. (a) FIG. 2. (a) Regions for the half- currents between 7 and the separatrix. (b) Decomposition of a self-intersecting circuit into the simple circuits Ti and F2. Open arrows show the cir- cuits orientations and simple ones the poloidal field direction. To study the properties of the equilibrium solutions, like the size or number of islands, we put n in terms of geometric parameters. For this, we develop a local successive approximations method that account for the non-nested topologies while keeping a simplified physical picture of the equilibrium. From ((T|) and Vx£ = fi j it can be verified that RV ■ {R- 2 Vip) -MOJ0- (5) Also, from the single-fluid MHD equilibrium [ll| the force balance condition Vp = j x B, leads to the relation IMiU = MoiVWO + R^Fi^F'Oj), (6) where the prime denotes d/dip, p(ip) is the kinetic pressure and F(ip) is defined after (pj. Equating (0 and ((6]) give s the well known Grad-Shafranov (G-S) equation 12|, |l3(, a nonlinear elliptic partial differen- tial equation. In ©, self-consistent surface functions {p(ip) , F {ip)} and boundary conditions characterize the magnetohydrodynamic equilibrium reached by the sys- tem 11] . In regular situations the equilibrium is de- scribed by a single family of magnetic surfaces, requir- ing a single choice of {p(?p),F (i/OI an d a single bound- ary (the plasma edge). This may not be appropriate for CRECs, but in a first approach a single choice of the sources p(tp) and F(ij)) can present current density inver- sions and multiple magnetic families [5|, l6j . Other works uses successive approximations to the solution with pre- scribed zero-order current density models and boundary conditions @,ll]. Formally, a broad range of choices for the arbitrary functions p(tp) and F(ip) or zero-order profiles may lead to current reversals and non-nested configurations, but the underlying description of the equilibrium topology is not restricted to the particularly chosen model. In the following we assume directly the existence of a small negative minimum of the current density in some small region of the plasma. With this approach we cover a wide range of arbitrary choices without specifying the global configuration of the plasma. Near a current density minimum at ro, a usual Taylor expansion gives j^r) « j^ro) + l/2(8r ■ V)l„j<p(r) + 0(Sr 3 ). For an up-down symmetric equilibrium this can be casted like j^ = jo + u{r 2 — nr 2 cos 2 9), where (r, 6) are local polar coordinates centered at ro. Here, jo — j(r ) < 0, r = |i — r | and 9 is measured counterclockwise between en and r — Tq. The parameters i and k are related to the curvature and ellipticity of j(r) about r . For k — 0, a circle with radius a = \f— 2j /l contains a zero current and for r < a/y/2 we have j^ < 0. Since we are interested in the topology of the magnetic surfaces in a region containing an almost vanishing toroidal current we define a region of interest with length 2a, centered at ro (Fig®. Q region of interest ( ) plasma I U < o FIG. 3. Representation of the region of interest. Provided that the region of interest is small we can solve ([5]) by a local scheme analogous to the successive approximations method 11]. In our case the inverse as- pect ratio, e = a/Ro, is a reliably small parameter defined from the region of interest instead of the plasma radius, in consequence, a fast convergence of the approximations is expected. Since the local current description is not in- tended to provide the vanishing of the current density at the plasma edge, the following solution is valid for r < a. Defining R = R + ax and z about the current minimum like <> d y-T ex ay we can write (O 9 I U = (l+et)(l-2r i +2/t4 (7) with r measured in units of a and ip in units of fj,o\jo\a 2 Ro- The parameter k corresponds to the squared eccentricity of the current density level sets. The solutions of ([7} will depend on two parameters only, allowing a detailed study of the bifurcations that change the topology of the equilibrium in the region of interest. Assuming a large aspect ratio e << 1 and small el- lipticity |k| << 1, we write the nondimensional flux as il){r,0) = 4>o(r) + eipi(r,x) + 0(e 2 ). Excluding the el- lipticity in the zero-order calculations we have ipo( r ) — (1 — r 2 /2)r 2 /4, consequently, the first order problem be- comes AV>i 3 - hr 2 2 k (8) with A the ing Tpi(r,x) Laplace operator for a plane. Dcfin- = tpi(r,x) + nx^/Ge we have Ai/ji = (3 — 5r 2 ) x/2. Then, we reduce the problem to a sim- ple ODE by introducing the ansatz ipi = fi r )x- After straightforward integrations and substitutions we obtain the poloidal flux function to a first order in e. rb{r,x) = 1 - -r 2 ) r — + ex [ 1 - % 2 16" + 6* ■ (9) With this local solution we can study the magnetic field in the region of interest r < 0(1). For instance, if wc write the toroidal field as B t — Bq/(1 + ex), the mag- netic pitch angle, defined by 7 m = arctan(_B p /_B t ), be- comes 7 m = arctan(r|V?/'|), with r = /io a |Jo|/-E>o- By a suitable choice of the region of interest and adjust of k we can fit experimental measurements of j m about a current minimum and model the corresponding magnetic surfaces (Fig. [4]). This local solution agrees well with the 0.023: (cm) : 3.4 R{m) -9 -6-30369 FIG. 4. Reproduction of experimental magnetic pitch angles in a current hole and reconstruction of the magnetic surfaces. Experimental data is from N.C. Hawkees et at, Phys. Rev. Lett. 87, 115001 (2001). reported Motional Stark Effect measurements in a cur- rent hole in [3j and predicts jo ~ —42.6 kA/m 2 , a value within the error reported by Hawkees et al. The structure of the magnetic surfaces depend on the relative positions of the critical points (VV> = 0), their type (elliptic or hyperbolic) and the value of ip on them. Since r and x are independent coordinates for ^ {0, 7r}, the off-axis critical points satisfy d r ip = d x ip = 0. The origin r — 0, is always an elliptic point between two crit- ical points in the x-axis with x c ~ ±1. The type of the these points depend on the ellipticity K. To under- stand this, we use the equation resulting of d r ip — and the condition for off-axis critical points \x c \ < r c . Tak- ing into account the domain of the first term in (0) we obtain |k| > e/8. The increase of the ellipticity leads to bifurcations about |«| ~ e/8, where off-axis critical points are created from those in the x-axis that change of type after the bifurcation occurs. A negative (positive) ellipticity represents a vertically (horizontally) elongated current density. For k < e/8, there is a single positive current channel about the negative one (Fig. [SJi-c). For K < ~ e /8 the positive channel contains an internal sep- aratrix (Fig. [5k,b) that for large anisotropies (or small e) k << —e/8 tend to split the channel by merging with the main separatrix (Fig. (5^,). For —e/8 < n < e/8 no off-axis critical points exists (Fig. |5fc) and for k > e/8, the on- axis hyperbolic point changes to an elliptic one creating a second positive channel (Fig. [5ji). -0.3e- -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 r\ does not correspond to a vanishing ellipticity. This is due to the toroidicity (represented by e), that adds an implicit anisotropy to the system. For this local solu- tion the anisotropy only vanishes for the cylindrical case (e = 0, k = 0) and never becomes negative, so it verifies i+ > 21- . The value of r\ saturates during the bifurca- tion at k — e/8 where, briefly, a zero current island is created. If we consider a very small region of interest (e — > 0), the current ellipticity dominates and the equi- librium always presents two islands. This means that, if a quasi-stationary transition to an internal reversed current density is possible, the magnetic axis is simulta- neously split into three separated axes. FIG. 5. Levels sets of i[> from Q with e = 0.1. The continuous separatrix defines the current channels and the dashed one the island internal structure. FIG. 6. Change in the anisotropy r\ as a function of the el- lipticity k for different values of e (left). Intrinsic anisotropy due to the toroidicity with zero ellipticity (right). In summary, the magnetic topology related to current density reversals defines several current channels within the plasma. The ratio between the current in the positive channels and the central negative current depends on a topological parameter measuring the anisotropy of the positive channels. In general terms the positive current is about twice the size of the central negative current. This anisotropy was shown to be related to the geometrical properties of the equilibrium in a region of interest inside the plasma. Then, it was observed that the anisotropy is always positive, indicating that the positive channels have at least twice the current magnitude of the negative one, causing the screening of this channel and forming a structure with net positive current. Finally, it was also shown that experimental magnetic pitch measures in a current hole are compatible with the existence of a small negative value of the current density. This work was partially supported by Conselho Na- tional de Desenvolvimento Cientifico e Tecnologico and Fundagao de Amparo a Pesquisa do Estado de Sao Paulo. Since the ellipticity k, and inverse aspect ratio e, de- fine the island sizes and the equilibrium topology, they must be related to the anisotropy r], in a continuous way. This is depicted in Fig. [6] Increasing the vertical elon- gation of the current density leads to a growth of the anisotropy of the current channels. The same is not true for the horizontal elongation, since the minimum * davidcirotaborda@gmail.com t ibere@if.usp.br [1] B. J. D. Tubbing, N. A. C. Gottardi, B. J. Green, J. A. How, M. Huart, R. Konig, C. G. Lowry, P. J. Lomas, P. Noll, J. J. O'Rourke, P. H. Rebut, D. Stork, A. Tanga, A. Taroni, and D. J. Ward, Nucl. Fusion 32, 967 (1992) [2] J. Li, J. Luo, S. Wang, P. Fu, B. Shen, F. Liu, B. Wan, [6] J. Shan, G. Xu, J. H. J. Yu, J. Hu, Q. Yuan, and Y. Hu, [7] Nucl. Fusion 47, 1071 (2007) [3] N. C. Hawkes, B. C. Stratton, T. Tala, C. D. Chal- [8] lis, G. Conway, R. DeAngelis, C. Giroud, J. Hobirk, E. Joffrin, P. Lomas, P. Lotte, J. Mailloux, D. Mazon, [9] E. Rachlew, S. Reyes-Cortes, E. Solano, and K.-D. Zas- trow, Phys. Rev. Lett. 87, 115001 (2001) [10] [4] T. Fujita, T. Suzuki, T. Oikawa, A. Isayama, T. Hatae, O. Naito, Y. Sakamoto, N. Hayashi, K. Hamamatsu, [11] S. Ide, , and H. Takenaga, Phys. Rev. Lett. 95, 075001 [12] (2005) [13] [5] A. A. Martynov, S. Medvedev, and L. Villard, Phys. Rev. Lett. 91, 085004 (2003) S. Wang, Phys. Rev. Lett. 93, 155007 (2004) P. Rodriguez and J. Bizarro, Phys. Rev. Lett. 95, 015001 (2005) C. G. L. Martins, M. Roberto, I. L. Caldas, and F. L. Braga, Phys. Plasmas 18, 082508 (2011) The separatrix must be followed so that the trajectory is differentiable This may be generalized to include several elliptic points and some internal separatrix J. P. Freidberg, Rev. Mod. Phys. 54, 801 (1982) V. D. Shafranov, Sov. Phys. -JETP 26, 682 (1960) H. Grad and H. Rubin, in Proceedings of the Second United Nations International Conference on the Peace- ful Uses of Atomic Energy, Vol. 31 (Geneva, 1958) p. 190